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1 INTRODUCTION 



Scalar-tensor theories (|Brans et ail ll96ll ) are among the the- 
ories able to explain the Universe late time acc elerated 
expansion (|Riess et all ll998f)(|Perlmutter et a] I1999T) . They 
generalise the ACDM or quintessence models by assuming 
a varying gravitation function G. Consequently they can 
be difficult to constrain because one has to interpret ob- 
servations, in particular cosmological ones, by consi dering 
the v ariation of G. Among the various experiments ljChibal 
l201lt> attempting t o constraint G, lunar laser ranging 
I Williams etaJ 2004) gives some stringent limits on the vari- 
ation of dG/dtG" 1 today. Another constraint comes from 
supernovae data whose magnitude-redshift relation have to 
take into account the variatio n of G on the supernovae 
mechanism (|Riazuelo et a 3 120021 ). 

In this paper we use lunar laser ranging and supernovae data 
at low redshift (z < 0.1) to constrain scalar-tensor theories. 
We show that these two kinds of observations are comple- 
mentary. Low redshift data allow to consider power series 
of the field equations near z — 0. This leads to constrain 
observationally some of the terms of these series without 
specifying the form of the scalar field potential U. It is then 
also possible to derive some constraints on the free parame- 
ters of a scalar-tensor theory once U is given. In particular, 
we will choose to consider the theories defined by G _1 = <j> 
with U = A, l/2m 2 (f) 2 and me A< *, A and m being the free 
constant parameters. 

The plan of the paper is as follows. In the first section, we 
write the scalar-tensor theory field equations in two ways. 
The first one depends on G and its derivatives and allows 
to use the lunar laser ranging. The second one consists in 
rewriting the field equations as General Relativity with a 
perfect fluid modeling the scalar field and defined by a den- 
sity p, a pressure p and an equation of state w. It is more 
appropriate to use the supernovae data. Developing some 
functions (scalar field, Hubble function, etc) as power series 
of the redshift near z — 0, we express some of their constant 
coefficients with observational quantities. In the second sec- 
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tion, we use lunar laser ranging and supernovae data to con- 
strain these coefficients. In a third section, these results are 
discussed and applied to some classes of scalar-tensor theo- 
ries. 



2 FIELD EQUATIONS 

The scalar-tensor theory action writes 



S = 



U + 16nL„ 



<j> is the scalar field, G(4>) the gravitation function, oj((j)) the 
Brans-Dicke function and U(<j>) the scalar field potential. L m 
is the matter Lagrangian. We will choose units such that 
G(to) = Go = 1 with to the proper time today. Varying the 
action with respect to the metric function g M „, we get the 
field equations for a flat Universe: 
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We want to rewrite them as General Relativity with a 



perfect fluid (|Capozziel"iol 120061 ) . 

above equations as 

1 



For this, we transform the 
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and we define the impulsion-energy tensor of the perfect 
fluid as 

= -(G -1 - l)(Ru.v - -^gnuR) + ^0, M </v - 
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We checked that T^ v ^ is conserved by calculating T?K. . 
We then recover the Klein-Gordon equation that is usually 
obtained by varying the action with respect to 0. We define 
the density p and pressure p associated to T^u,) as 

P = T WW = -HG- 1 - 1)H 2 + \^4?- 3H(G~ 1 )- + \U{1) 



H (2p - Pl -U(l)-3U^(1)) 



6# 2 (4p 



2U(1) - 



31^(1)) - (2p 



(4) 



with (7(1) = [7(0 = 1), the present value of the potential, 
U<f,(l) = ^{<t> = 1) and Q = 36Hi + 4p - 2U(1). In the 
next section, we constrain 0i an d wo thanks to lunar laser 
ranging and low redshift supernovae and then Hi, U(l) and 

02- 



P = T uw /a 2 = (G- 1 -1)(2H + 3H 2 ) + ^ 2 + 

(g- 1 )- + 2H(G- 1 Y - hj 



where a dot means a derivative with respect to the proper 
time t. Hence, the equation of state w = p/p is 



1 



-H(G~ 



21G' 1 -1)H + (G" 1 )" + f0 2 



-3(G-! - l)W + if 2 - SHiG- 1 )- + \V 



This last expression agrees with the conservation of T^ v ^ 
since it checks p — —3Hp(l + w). 

In the rest of the paper, we redefine the scalar field 
without loss of generality such as the kinetic term % 0, M 0' M 

is cast into 2 < P,^4 > '^ ■ We choose to study the class of theories 
defined by G _1 = <f>, leaving U unspecified. 

Around z = 0, we write the density p and the equa- 
tion of state w as some power series p = po + p\z + 0(2) 
and w — wo + w\Z + 0{2). Then, introducing these series 
in the energy conservation equation for p, we get for the 
zeroth order term 

pi = 3(1 + wo)po 

It is possible to get equations for any higher order terms but 
they introduce some constant coefficients p n and w n that 
we cannot constrain observationally. So we do not consider 
them here. Using this expression for p\ in the constraint 
equation 



H 2 = Hi 



n m0 (l + z) 3 + (l-fi m0 ) — 
po 



with 1 — Q.mO = Oo = and writing the Hubble function 

o 

as H = Ho + H\z + 0(2), we get for the first order term 
(linear in z, the zeroth order term giving the usual constraint 
f2 m o + fio = 1) 



H x 



-H [-3 + 3w (n m0 -1)] 



(2) 



Developing cj> as <f> — <f>o + <f>iz + (j> 2 z 2 + 0(3), it comes cf>o = 1 
since we choose Go = 1. Moreover, from the definition (fTJ) 
for the density p, we get for <f>\ and (j>2 considering zeroth 
and first order terms 



-6T 



Ho 



(3) 



3 CONSTRAINTS ON SCALAR-TENSOR 
THEORIES WITH LUNAR LASER 
REDSHIFT AND LOW REDSHIFT 
SUPERNOVAE 

In the rest of the paper we adopt the value Q m o = 0.27. 
The results of this section weakly depend on the value of 
fio but on the value of Ho- For this reason, we will choose 
for Ho the WMAP recommended value Ho = 71 ± 2.5. The 
value of 0i can be deduc ed from lunar laser ranging. Indeed, 
from lunar laser ranging (|Chiball201ll ). we know that today 
(G/G) t0 = (4 ±9) x HT 13 ^- 1 . But 

G/G = -(G" 1 )' /G" 1 = (G" 1 )'G^ 1 H(1 + z) 

with a prime meaning a derivative with respect to 
z. Developing G _1 = <j> in power series, we thus 
find (G/G)t = <piH in z = 0. It follows that 
4>i — 0.0055 ± 0.012. We will consider the plus sign 
for 0i in (J3J, and consequently the minus sign for 02 in @, 
corresponding to the positive best fit value of 0i got with 
lunar laser ranging. 

In what follows, we are going to use these values of 
Ho and cj>i as priors to determine Ho, (j>i and wo with 
supernovae data. For that, since G evolves with time, we 
have to modify the magnitude-red shift relation usually used 
in General Relativity. Following (|Riazuelo et all l2002h , we 
write 



m — 51nd; 



4 



with m th the (theoretical) modulus distance and d; = 
c(l + z) J Q Z H(z)~ 1 dz the usual luminosity-distance in a flat 
Universe. At small redshift, we get 



m th = (51nz— — In 0o+5 In -2— ) + (5- 



15 



5 Hi 



«+0(2) 



15 0i 

4 "~ r Ho' ' v " 4 0o 2# 

Since 0o = 1 and considering the form ((2} for Hi, this 
rewrites as 



c 15 5 

= (5 In 2+5 In— )+(5 — t-0i+- [-3 + 3w (Ono 



Ho' v " 4 

Let us remark that 0i and u>o both appear in the linear 
term. It means that supernovae data can only constrain a 
combination of these parameters and not each of them sepa- 
rately. The bound on 0i got with lunar laser ranging is thus 
necessary to constrain efficiently wo with supernovae data. 
To proceed, we use the standard x 2 minimisation with some 
priors on 0i and Ho, i.e. 



l)])z+0(2) 



216ijp 4 ± y/g [36ff 3 + Hi{2 Po - [/(!))] 



+ 



2 

X = 



E 

P =i 



(m° bs - m th ) 2 



+ 



0.0055) 2 (Ho - 71) 2 



0.012 2 



2.5 2 
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n is the number of supernovae and m obs their observed dis- 
tance modulus. We consid er the 166 supern ovae with z ^ 0.1 
from the last Union data (|Amanulrahll2010l ). 
At la, one finds that H = 69.4 ± 1.2, fa = 0.0055 ± 0.0187 
and wo = —0.53 ±0.68. We got similar results for fa and wo 
when marginalising Ho instead of considering it as a prior. 
Let us also remark that the A.CDM model defined by fa — 
and wo — — 1 is in agreement with the data. From these val- 
ues, we also derive from Q that 

(7(1) = (2.20±0.06)10~ 35 s~ 2 

This value of the potential today (U(l) = U(fa)) is of the 
same order as the cosmological constant. From © we then 
deduce that 

Hi = (2.06 ±1.67) x 10~ 18 s _1 
and finally from @ that 

fa = 0.50 - 9.04 x lO 31 t/ (l)± 

(0.54 ± 1.15 x 10 31 (/^(1) ± 9.43 x 1O 64 ?7 (1) 2 ± 
8.17 x 10 63 a^ (1) ) 1/2 S " 2 

with 0u (15 the error on U^(l). Note that all the values 
derived in this section are independent on the form of U (fa 
but fa. In the next section, we discuss about these results. 



4 DISCUSSION 

In this paper we study a class of scalar-tensor theories de- 
fined by G~ x = 4> an d U — U(fa. Considering low redshift 
data allows then to get some constraints on the scalar field 
(i.e. wo, fa, U(l) and a constraint on the variation of the 
Hubble function with Hi) independently on the form of its 
potential U(fa by using lunar laser ranging and supernovae 
data. These constraints for this Brans-Dicke-like scalar field 
considered as a dark energy thus come from cosmology, at 
scales much larger than solar system scales. 
We note that low redshift supernovae tend to lower the value 
of Ho recommended by WMAP although it stays in the 
la confidence contour of this last experiment. The positive 
value of Hi shows that the Hubble function is decreasing 
today. Lunar laser ranging determine the value of fa and 
this allows to constrain wo with the supernovae data. With- 
out lunar laser ranging, supernovae data only constrain a 
combination of fa and wo- Constraints on u>o are rather 
weak with respect to what we usually get when one consid- 
ers high redshift supernovae but future project like JDEM 
should improve this situation. The constraints we got on wo 
and fa are in agreement with the ACDM model (for which 
wo = —1 and fa — 0). (7(1) is thus accordingly close to the 
cosmological constant value of the standard ACDM model. 
Constraints on fa depend on the form of the scalar field po- 
tential. 

Let us apply the above results to some specific scalar-tensor 
theories. The direct generalisation of the A CDM model for a 
scalar-tensor theory corresponds to U = A( Mas lankall 1983^ ) . 
For this model, we get A = 17(1) = (2.20 ± 0.06)10 _35 s- 2 
and fa = 0.50 ± 0.74. This last range of values of fa is in 
agreement with the standard ACDM model since it contains 
the value fa = 0. Another interesting model is based on a 



varying potenti al U = I /2m 2 fa where m can be interpreted 
as a mass term (|Lindell2005l ). The constraint on U(l) allows 
to deduce that m 2 = 2U(1) = (4.4 ± 0.12)10~ 35 s -2 and 
then fa = 0.50 ± 0.74. This value of fa is similar to the 
previous model, showing its robustness despite the possible 
variation of U. It could change if we consider a mo del with 
two free parameters like U = me A ^ (|Halliwelll ll986T ) . Then, 
we get me A = U{\) = (2.20±0.06)10~ 35 s~ 2 and fa = 0.50- 
0.0019A ± V0.54±0.25 x 10" 3 A±0.45 x 10" 4 A 2 . Here we 
can only constrain a combination of parameters i.e. A and 
m. But if A >> 10 2 , fa has still the same value as with the 
two previous models and is still in agreement with a ACDM 
model. 

The method of this paper can be applied to other scalar- 
tensor theories with various forms of G. A larger number 
of low redshift supernovae and a better determination of 
Ho should allow to improve the constraints we got and to 
exclude or not a variation of the scalar field via the deter- 
mination of wo, fa and (f>2. 
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